Magnetism in Closed-shell Quantum Dots: Emergence of Magnetic Bipolarons 
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Similar to atoms and nuclei, semiconductor quantum dots exhibit formation of shells. Predictions of magnetic 
behavior of the dots are often based on the shell occupancies. Thus, closed-shell quantum dots are assumed to 
be inherently nonmagnetic. Here, we propose a possibility of magnetism in such dots doped with magnetic 
impurities. On the example of the system of two interacting fermions, the simplest embodiment of the closed- 
shell structure, we demonstrate the emergence of a novel broken- symmetry ground state that is neither spin- 
singlet nor spin-triplet. We propose experimental tests of our predictions and the magnetic-dot structures to 
perform them. 
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Formation of shell structure is a ubiquitous feature in fi- 
nite fermionic systems, such as atoms, nuclei, and quantum 
dots (QDs) [1]. An effective potential, in which the fermions 
are assumed to move independently, can be attributed to the 
underlying mean field, arising from an interplay of particle- 
particle interaction and confinement. Open- shell atoms, e.g., 
Ag and Fe, undergo a spontaneous symmetry breaking of 
the mean field and are magnetically active due to their spin- 
polarized ground state (GS). For open-shell QDs doped with 
transition-metal atoms, typically Mn, strong exchange cou- 
pling between a carrier spin and the impurity spin is expected 
||2|-l5|], Such QDs exhibit magnetic ordering, which persists 
even up to room temperature ||6|4sl]- In contrast, closed- shell 
fermionic systems, e.g., noble gases, are known for their sta- 
bility and the total spin-zero GS, making them magnetically 
inert §]. 

According to a theorem by Wigner [10], the GS of any non- 
magnetic two-electron system, including a two-electron QD, 
is a spin-singlet. Thus, it would seem that closed-shell QDs 
doped with Mn do not allow magnetic ordering. However, on 
the example of a two-particle (two electrons or holes) system, 
we show that the Mn doping does alter the magnetic properties 
of closed-shell QDs. Surprisingly, we find a GS, which is 
neither a singlet nor a triplet, and allows ordering of Mn-spin, 
owing to the spontaneously broken time-reversal symmetry 
[11]. This mechanism of magnetism is different than in the 
open-shell systems, such as bulk (Ga,Mn)As or (Cd,Mn)Te 
11121 . By definition, the open-shell systems have more of either 
"spin-up" or "spin-down" carriers. This is in contrast to the 
magnetic closed-shell state considered here, characterized by 
zero total spin projection. 

Carriers confined in a QD interact with magnetic ions via 
contact exchange interaction, described by 



H ex = -(Jex/No) £ s ia g a pS jpS(ri-Rj), 
ij,ap 



(1) 



where No is the cation density. The exchange integral, J ex , 
is typically ~ 0.1 eV for electrons, and ~ — 1 eV for holes. 
Carrier and magnetic ion positions are denoted by r and ir- 
respectively; S and s are the Mn and carrier spins, g-tensor 



describes possible exchange-coupling anisotropy, which is 
caused by spin-orbit interaction combined with the quasi-two- 
dimensional shape of QDs. In many semiconductors, this 
anisotropy is almost negligible for electrons. In contrast, 
for confined holes, the spin-orbit coupling leads to a strong 
anisotropy with "easy axis" along the growth direction z 111311 . 
Thus, Eq. (OQ) reduces to the Ising Hamiltonian for the heavy 
hole ±3/2 pseudospin subspace. 

We focus on a two-carrier QD, the simplest example 
of a closed-shell system. The total Hamiltonian, H = 
Hf ± H^, contains the fermionic part Hf, which employs 
a typical two-dimensional (2D) model for two carriers in a 
QD 0, H f = -h 2 / (2m*) (Vf ± V*) + nf(0 2 (r 2 + r 2 ) /2 + 
e 2 1 {Aner\ 2 ), where h is the Planck constant, m* the effective 
mass, C0o determines the 2D confinement, e is the electron 
charge, £ the QD's dielectric constant, and ryi = \r\ — r 2 \. 
The Coulomb interaction is characterized by the effective Ry- 
dberg energy Ry* — nfe A /[?>2 (nek) 2 }. We express the two- 
fermion wavefunction <I> as 



(2) 



Here XaiX) an d XoQ) are spinors of the carriers 1 and 2, 
and (7 = ±1/2 (or t» 40 correspond to the spin projection 
along the quantization axis z. The Pauli principle requires 
<Poo'(n,r 2 ) = -<Pc7'cj02,n). Thus, (pft and (p u must be 
antisymmetric functions of r\ and r 2 , while 1,7*2) and 
<P;t( r 2? r i) transform into each other and could be neither 
symmetric nor antisymmetric. To understand the origin of 
these states, we first consider the "non-magnetic" Hf, spin- 
independent and invariant under the r\ <-> r 2 interchange. 
Thus, its eigenstates are also eigenstates of the total carrier 
spin, and (p GG i {r\,r 2 ) are either symmetric or antisymmet- 
ric. The GS of any system of two identical spin- 1/2 fermions, 
described by a spin-independent Hamiltonian, is a singlet 
with c + o' = and (p^i(r u r 2 ) = (p^(r 2 ,ri) |Tfl[l5|]. Ex- 
cited states are either singlets, or triplets with (p a(7 f (r 1,1*2) = 
— (p G(J f (r 2 , r\). Symmetry of the singlet (ground- state) and 
triplet states is illustrated in Fig.Q] 

The above classification does not fully apply when H ex is 
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FIG. 1 . (color online) Two-fermion states described by the symme- 
try of the orbital and spin part of their wavef unctions. Arrows show 
spin projections (up or down), while radii of the spheres indicate the 
extent of orbitals. For the S state, the radii of the orbitals correspond- 
ing to spin-up/-down are the same. For PS, the spin-down orbital is 
larger, [Eq. (|3)], leading to finite spin density and a magnetically- 
active state. 




FIG. 2. (color online) Emergence of magnetic bipolarons in QDs. 
(a) Without carriers, Mn spins (light arrows) are randomly oriented. 
Double arrows in (b, c) show carrier's spin projection associated with 
the orbitals (solid and dashed lines), (b) With 1 carrier, a magnetic 
polaron forms, lowering the total energy of the system, due to the 
coupling of the carrier's spin with the induced Mn magnetization. 
Mn spins align in one direction, (c) Two carriers assemble in a PS 
state, forming what we term a magnetic bipolaron. The sign of the 
Mn-spin projection depends on the sign of the carrier- spin density 
(difference between dashed and solid curves). The extent of the or- 
bitals (length of dotted lines) is different. 



included. We first consider a two-hole system, assuming the 
Ising exchange, which allows to express any eigenfunction of 
H as \j/ = ({r/,5 1 /} ; {Sj z }) YljX J (Sjz)> where s t are the spin 
variables, Sj z is the spin projection of the j-th impurity, and 

' :X J (S jz )=Sj zX J (Sj z ), 
j Z = —J, . . . ,J with J = 5/2. This separation of \jf (or 



X J (Sjz) is an eigenfunction of Sj z 
where S 

"classical approximation" for Heisenberg spins) is correct to 
order N^/ 2 , where A^in is the number of Mn spins fl3l . It 



has been widely used in the literature H6HJ/ZD- The separation 
is exact in the Ising case, so that approximations, such as the 
variational method used below, are needed only for the carrier 
subspace. The two-hole depends on {Sj z } parametrically, 
and it is this coupling that leads to formation of magnetically 
ordered states. The z-projection of the total carrier spin, Z, 
is a good quantum number, so that the Hilbert space splits 
into three orthogonal subspaces with values £ = a + a' = 0, 1 



and —1 lH3l . We show that the GS is never a singlet, i.e., 
9(7(7' (^1^2) 7^ 9(7(7' (r2,r\). Instead, it is either a triplet (T) 
with |E| = 1, or what we term a pseudo-singlet (PS), with Z = 
0, which reflects its closed-shell character. Unlike the typical 
singlet, PS leads to ordering of the magnetic moments of the 
open-shell d-orbitals of Mn, due to breaking of time-reversal 
symmetry. 

Magnetic polarons form by aligning a "cloud" of Mn- spins 
by a single carrier localized in, e.g., a QD or in an impurity 
potential (l6lll9|]- The consequences of presence of two car- 
riers in a magnetic QD are shown in Fig. O We predict that, 
even when the PS is the GS, a magnetic bipolaron is formed, 
despite vanishing L. The finite exchange interaction is possi- 
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FIG. 3. Spin corral. Colored surface: The hole-spin density pps (arb. 
units) of PS. Black circle indicates pps (r = Rq) =0. Green arrows: 
Mn spins, placed at a radius Rq, which maximizes the stability of the 
ferromagnetic alignment. Red and blue arrows: The more probable 
hole- spin projections at two positions. The parameters: h(Oo = Aq = 
30 meV, Ry* = HcOq/IO, m* = 0.5m , w = 1 nm Q. 



ble due to different Bohr radii of the "up" and "down" orbitals 
of PS, see the lowest row of Figs. [T] and Oc). The non-zero 
exchange and symmetry breaking is particularly obvious for 
Mn spins arranged in a ring ("spin corral" ), Fig. [3l These 
spins are "ferromagnetically" ordered, their common direc- 
tion marks one of the two possible states (Mn spins point- 
ing either up or down, perpendicular to the QD plane), cor- 
responding to two stable magnetic-bipolaron solutions sepa- 
rated by an anisotropy barrier. The latter is defined by the 
strongly anisotropic hole g-tensor. 

To analyze magnetic-bipolaron states and the symmetry 
breaking induced by H ex , we approximate the GS of two in- 
teracting holes using two alternative trial wavef unctions. The 
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first one is the PS (£ = 0) 



1 



-u(r 2 )d(r l )xt(2)^Xi^)], 0) 

where d = ^2/7lL~ l d exp (^—^/L 2 ^ are single-carrier or- 
bitals corresponding to spin "up"and "down" respectively, 
and L u ^ are the variational parameters. Comparing this with 
Eq. ©, we find (p n (r u r 2 ) = 2- lt2 u(r x )d(r 2 ) & The 
second is a Z = 1 triplet 



<D T = (p n (n,r 2 )^ t (l) <g>£ t (2), 



(4) 



2/(7rL: 



co- 



where (p n (r u r 2 ) = -(p n (r 2 ,r { ) 
[ne^ -r 2 e^ 2 ] eh( r i+ r 2)/4] i s the orbital part, in 
ordinates r, 0, and with one variational parameter Lj. 

The exact treatment of the isotropic g-tensor (for two- 
electron QDs) requires a large dimension of the Hilbert space. 
The problem can be circumvented by replacing the Mn-spin 
operators with classical spin vectors Biol . Then, the two 
PS and two T solutions (Mn spins pointing up or down, see 
above), found for holes, form continua (one for PS and one 
for T), corresponding to Mn spins aligned along arbitrary di- 
rections. Any particular solution preserves the form given by 
Eq. © [with the z-axis parallel to spontaneous magnetiza- 
tion], and it corresponds to formation of the magnetic bipo- 
laron. The distinct feature of the isotropic case is lack of an 
anisotropy barrier between different solutions belonging to the 
same continuum. 

Owing to the disk-like shape of typical self-assembled 
and vertical QDs, the z-dependent Schrodinger equation 
is factorized out, while the height, h, (along z) of such 
QDs is usually small. This allows to assume that only 
the lowest level of this equation is relevant B22I1 . Thus, 
for the PS state [recall Eq. ©], the matrix element of 
H ex is E ex = -J ex /(hN )Zj PPs(Rj)Sjz, where pj>s(Rj) = 
| u (Rj) | 2 — \d (Rj) | 2 /2. If Luj^Lj, then E ex ^ in a mag- 
netic QD. 

In general: For H ex = (non-magnetic QDs), PS reduces 
to a singlet for any Ry*, because the non-magnetic total en 
ergy functional E$ reaches a minimum, E$, for L® 
(0 indicates a variational minimum). In all studied systems 
with H ex ^ 0, however, u° ^ d°, i.e., PS does not reduce to a 
singlet. Because of its larger spin density, T may become the 
GS, despite its non-magnetic energy being higher than E$ . (T 
as the magnetic GS was discussed by Govorov 0,[23|]). The 
functionals £"ps,t reach minima, when Mn spins are antiparal- 
lel to hole- spin density. 

Because of the exchange coupling, Eq. (OD), the total energy 
functional contains a term linear in the carrier- spin density 
s(r). This term leads to instability of the closed-shell singlet 
with s(r) = 0. PS, on the other hand, must satisfy a weaker, 
integral constraint J s(r)dr = 0. Thus, any small variation 
of the wavefunction <I> that promotes a non-zero s(r) Pps, 



while preserving the integral constraint, leads to the instability 
of the singlet state, because the variations of the kinetic and 
the Coulomb energies contain only second or higher powers 
of s(r). 

We now describe our results for two particular distributions 
of Mn spins, starting with homogeneous Mn content x. Fig- 
ure a) shows the phase diagram of PS and T. The former re- 
mains the GS for moderate values of Ry* and of the saturated 
exchange splitting Aq = Jx\J ex \ [24]. We analyze the GS by 
considering small variations, 5, of the characteristic lengths 
from their non-magnetic values L<! and Lj nm . We write 
L u ,d = L° s (1 ± 5) for PS, L T = L° nm (1 + 5) for T, and treat 
8 as a variational parameter, see Figs.^b), (c), and (d). The 
quantity AEps (<5) = E?$ (8) — E$, plotted in Fig.[4tb), has two 
minima corresponding to the two opposite Mn-magnetization 
profiles mentioned above. The E^ s — E<! gap can be an order 
of magnitude larger for colloidal QDs with a few nm diame- 
ter B25I1 . suggesting stability of PS at liquid nitrogen temper- 
atures. Figure HJc) shows that AEj(S) = Ej(S) — E® has S 2 
dependence with a single minimum. Unlike PS, the T wave- 
function at the variational minimum is the same as for H ex = 
(i.e., 5° = 0), while its energy is lowered by Ao. Hole- spin 
densities and Mn-spin profiles corresponding to PS and T are 
shown in Figs. Hid) and (e). The small variation of <I> dis- 
cussed above is °< 8. The singlet-PS instability manifests itself 
as the cusp in the solid line in Fig.Htb). 

Results of a full variational calculation [20], [e.g. Fig.lHa)], 
confirm the validity of the 5 approximation. By studying 
the radius at which M z changes sign, we deduce the mag- 
netization profile for an inhomogeneous x distribution, such 
that x(r < Rq) = 0. Independently of other details of the in- 
homogeneous distribution, pps has the same sign at all Mn 
sites, leading to a "ferromagnetic" alignment in a large class 
of QDs. The same alignment arises when T is the GS. 
CdSe/(Zn,Mn)Se epitaxial QDs with Mn only at the periphery 
were created by intentionally introducing Mn in the material 
surrounding the dot |2ll26|]. The placing of individual Mn ions 
with a scanning tunnelling microscope B27I1 . is a promising 
path to realize the "spin corral". Additionally, x (r < Rq) = 
could be realized in colloidal QDs, where radial segregation 
of impurities occurs during growth |25|]. Such systems show 
strong exchange coupling |6j], and can be controllably charged 
|7|], thus avoiding fast Auger decay in type-I QDs with two 
electron-hole pairs. This decay can also be suppressed using 
core-shell colloidal nanocrystals, equivalent to type-II epitax- 
ial QDs, due to electrons-holes separation 1 81 12811. We estimate 
that for colloidal QDs with ~ 5 nm diameter 112511 . the singlet- 
triplet splitting can be ~ 100 meV, resulting in a PS ground 
state for a wide range of parameters. 

PS existence can also be experimentally verified for 
homogeneous-x QDs, e.g., as a blue shift of interband photo- 
luminescence with magnetic field, B, applied along the z— axis 
1 12]. With increasing B, all the Mn spins will tend to align 
antiparallel to it, destroying the PS magnetization profile, 
Fig.HJd). Thus, PS should increase its energy, while evolving 
towards the ordinary singlet. This effect could be observed 
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FIG. 4. (color online) (a) Phase diagram of two-hole ground 
states for homogeneous distribution of Mn. Horizontally (vertically) 
hatched area indicates the ranges of Ry* and Aq, for which the ground 
state is PS (T). PS reduces to singlet only for Aq = 0. (b) for PS, and 
(c) for T show the non-magnetic (dashed) and total (solid) energies 
as a function of 5. Panels (d) and (e) show the hole- spin densities 
(at the corresponding variational minima), along any direction in the 
x — y plane, with arrows indicating the Mn-spin profiles. The pa- 
rameters for (b)-(e) are Hcoq = 30 meV, A = \ho\), Ry* = hcoo/10, 
m* = 0.5m . Vertical line in (d): R = L° s /y/2. 



in type-II QDs, where the electrons (unlike the holes) reside 
in the barrier and do not modify the physical picture. Strong 
exchange coupling in these QDs is seen as magnetic-polaron 
formation |8|]. 

Our results can be generalized to closed-shell QDs with 
more carriers, and to systems, not described by Hamiltonian 
Hf, typically used for self-assembled 112911 . vertical 113011 . or 
lateral Q QDs. Shells also form for other confinements fll, 
e.g., colloidal QDs can be approximately described by spher- 
ical or ellipsoidal potential Bill . Closed- shell (L = 0) states 
appear even for asymmetrical QDs with many carriers ll32ll . 
More general choices of trial functions could further stabilize 
the PS energy, in some cases accompanied by lowering the 
symmetry of the Mn-spin alignment. 

This work was supported by DOE-BES, US ONR, AFOSR- 
DCT, NSF-ECCS, and CAREER. 
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